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Ch '. Abstract 

The best method of measuring CP violating effect in neutrino oscillation experiments is to 
construct and use a neutrino beam made of an ideal mixture of v e and v e of monochromatic 
lines. The conceptual design of such a beam is described, together with how to measure the 
£T) • CP-odd quantity. We propose to exploit an accelerated unstable hydrogen-like heavy ion in a 

storage ring, whose decay has both electron capture and bound beta decay with a comparable 
| fraction. 
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It is of great current interest, both to microphysics and to cosmology, to clarify how the matter- 
^ ■ antimatter asymmetry is generated. One key idea is the leptogenesis [1] in which the lepton asym- 
^ ■ metry of the universe is first generated by CP violating lepton number violating processes, and is 
Ph! later converted to the baryon asymmetry via high temperature B- and L-violating electroweak pro- 
r-| ! cesses. To advance this idea further, it is important to deepen our understanding on the origin of 
^ ; CP violation in the leptonic sector. 

Measurement of CP violation in the neutrino sector is thus one of the most fundamental problems 
facing physics beyond the standard particle physics. In this note we develop a concept of CP-even 
neutrino beam, which serves this purpose. Moreover, it is one of the most important experimental 
issues in neutrino physics to determine the relevant CP parameter in neutrino oscillation experiments. 

The ideal neutrino beam for CP measurement would be a mixture of monochromatic v e and v e 
beam for which the detector response is symmetrical, producing after the oscillation equal numbers 
of /i ± if CP is conserved. Such a beam is referred to as a CP-even neutrino beam in this paper. 
Our proposal is to simultaneously use the monochromatic neutrino v e of electron capture (EC) and 
anti-neutrino v e of bound beta {bf3) decay from hydrogen-like heavy ions. Such ions do exist, for 
example [2] 

114x48+ 110 a 46+ 108 a ^46+ 104 RV ,44+ /-, \ 

49 ln > 47 A S > 47 A S > 45 Kil ■ V 1 ) 

The concept of CP-even neutrino beam is realized, to a large extent, by this kind of beam. It is an 
extension of the ideas of monochromatic neutrino beam using EC proposed in [3], [1], and beta beam 
proposed in [5], [6]. 
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The hydrogen- like heavy ion of this kind has b(3 channels in addition to continuum beta (c/3) 
decay and EC, and thus produces monochromatic v e and v e beams once accelerated in a storage 
ring. One can arrange all lines and a part of the continuum neutrino energy to fall into an optimal 
range for detection: above HOMeV of the muon production, but below multi-pion production in 
detector placed at a distance. We shall argue that this provides an excellent opportunity of precision 
experiments to determine the CP violation phase S [7] as well as the mixing angle #13. 

Although not exhaustive by any means, a systematic method of search for candidate ions is briefly 
described, starting from data of neutral atoms. Since we momentarily ignore c/3 contribution for de- 
tection, one important measure for the CP-even beam is a large EC rate of neutral atom. This rate 
is almost the same for the case of the heavy ions. It is then important to look for calculated b/3 rate 
of about twice of EC rate, since cross sections for v and v has this ratio (for more details, see below). 

A possibility of the bound beta decay that produces monochromatic neutrino has been considered 
theoretically [8] , and their dramatic example that becomes possible only for highly ionized atoms has 
been demonstrated experimentally |9j. Relative strength of different neutrino lines is proportional 
to the atomic wave function squared at the nucleus. This factor for the s-wave state of the principal 
quantum number n is |^ ns (0)| 2 = (Z/nas) 3 /7r (using for simplicity the solution of the Schrodinger 
equation) for hydrogen-like atoms of charge Ze. The ratio of the bound to the continuum contribution 
is given by 

r B = ^N n (^^n^r^?K-\^), (2) 

with N n the multiplicity factor of available levels and Abe the difference of binding energy in neigh- 
boring ions. Here K(x) is related to the phase space integral of the continuum contribution. When 
the effect of Coulomb distortion of wave functions (the Fermi integral) is ignored, this function is in 
the range, K(x) = 0.15 ~ 0.74 for x = 1/2 ~ 2. When variation of the level difference Q c p + A B e 
with n is small, the fraction of Is contribution is 1/(2^(3) — 1) ~ 0.71, 2s contribution 0.18, and the 
rest 0.11. To obtain a large b/3 rate, r B = 0[1] is required; for instance, with Q c/3 + A B e = m e, the 
bound ratio r# = 1 corresponds to Z = 57. We thus look for heavy atoms of large Z. 

We found the best candidate ion to be 4° 8 Ag 46+ , and the next best one is 4g 4 In 48+ . The neutrino 
and anti- neutrino energy spectrum from decay of 4y 8 Ag 46+ at rest is shown in Figure 1. Effect of 
distorted plane wave under the nuclear Coulomb potential is important and included in the form of 
the Fermi integral. 

Table 1 gives their monochromatic energies and the endpoint energy of the continuum, along 
with their fractional contributions. To compare the line contribution with the differential continuum 
spectrum, we divided the line strength by an energy resolution factor AE, which is taken lOkeV in 
Figure 1. We denoted decay contributions to excited daughter nuclei by * in Figure 1 and Table 1. 
It is important to use the wave function of the Dirac equation rather than the Schrodinger equation, 
since the wave function at nucleus is sensitive to the short distance behavior. Most conspicuous lines 
from 108 Ag 46+ have fractions; 0.020 adding bound K-, L- and M-beta decays all to the ground nuclear 
level, and 0.0096 for EC + EC*. 

All promising candidates we know of are listed in Table 2, along with their characteristic features. 

CP violating effects are present only for the appearance experiment [10], hence we need to 
detect the muon neutrino from the oscillation, u e — ► and v e — > u^, thus to set the neutrino 
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Figure 1: Beta and EC neutrino spectrum from 108 Ag 



Table 1: Branching ratio of c/3, 6/3, and EC from ^Ag 464 " decay 



decay mode 


neutrino energy [MeV] 


branching ratio 


c(3 


1.649 


0.9479 


c/3* 


1.016 


0.0174 


6/3(ls) 


1.671 


0.0149 


bp(2s) 


1.648 


0.0041 


6/3(3s) 


1.645 


0.0012 


EC(ls) 


1.897 


0.0087 


EC*(ls) 


1.463 


0.0009 



Table 2: Neutrino energy (E) and branching ratio (BR) of cf3, b/3 and EC 



ion 


114 Tti 48+ 
49 111 


110 A ,,40+ 
47 A § 


half-life 


70.9[sec] 


24.4 [sec] 




E[MeV] BR 


E[Me\] BR 


cp 
6/3(1*) 
EC(ls) 


1.99 0.981 
2.01 0.012 
1.43 0.002 


2.89 0.947 
2.91 0.005 
0.87 0.001 




ion 


108 A„4(j+ 
47 A § 


1 4 u |Rh 44 + 


half-life 


2. 36 [min] 


42 [sec] 




E[MeV] BR 


E[Me\] BR 


c(3 
bp{ls) 
EC(ls) 


1.65 0.948 
1.67 0.015 
1.90 0.009 


2.44 0.971 
2.46 0.006 
1.12 0.002 




EC(ls) 
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energy above muon production. The appearance probability for v e — > (and v e 
approximation given by 



z/ M ) is to a good 



„2 -2 



9 9 9 SfTtloL 

) = S23 sin" 5 26>i3 sin"' 



4P 



+c 23 sin 26*12 sin' 



5m\ 2 L 



AE 



+Jcos ( ±8 



<5m 2 3 L\ 5m 2 L . Sm^L 



AE 



12 J 



AE 



sin 



n3 J 



4P 



(3) 



with J = C13 sin 26?i2 sin 26*23 sm 26* 
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Sm 2 2 /{Aml 3 ), which implies roughly #13 



We note that for the range of (<5m 2 2 /(4m 2 3 )) < 6^3 < 
i2/v*" t i37; wmu! impiico iuu 6 iu,y 1/13 = 0.0025 ~ 0.05, the last CP sensitive term in eq.Q is 
the largest. The matter effect that may mimic CP effects is negligible at this low energy. To de- 
termine #13 and 5 with precision, choice of the detector location L is important. Taking account 
of the neutrino flux factor oc 1/L 2 , three terms in the oscillation probability ([3]) have different 
L— dependence, oc L~ 2 , L° , L" 1 respectively, when the phase ip = 8m\ 3 L / (AE) is near the oscillation 
peak. To maximize simultaneously the flux at the detector and sensitivity to CP parameter 5, the 
best location is at the first peak ip = it/ 2. With this choice, the last term in (jSJ) becomes proportional 
to ±Jsin<5. Thus the symmetric combination of measured quantities P„ e „ +P PePfi is sensitive to #13, 
while the asymmetric combination P VeV — Pp e p M to the CP violation parameter J sin 5. This choice 
fixes the relation between E and L to be L/E ~ 517km/1000MeV. In Figure 2 we plot the oscillation 



probability P(V e 



as a function of 5 for three values of sin 2 26*13 = 0.1 , 0.05 , 0.01. 



0.08 



0.06 



> 
T 

o 

> 

Cm 



0.04 



0.02 




Figure 2: Oscillation probability P VeV at the first peak <p 
top to bottom). Pp c p M is obtained with 5 — > —5. 



tt/2 for sin 2 2^3 = 0.1,0.05,0.01 (from 



The neutrino energy E can be chosen at will by adjusting the acceleration 7. There are 
many practical factors to be considered; here we select them mainly from the viewpoint of neutrino 
detection. The cross section becomes larger as the energy goes up, but distinction between /i + and 
jjT , mandatory to establish genuine CP violating effects, becomes harder. Below the single-pion 
production threshold, backgrounds due to 7r ±7 s are absent; but the cross section is small and the 
inevitable Fermi motion smears out kinematic relation between /i^ energy and scattering angle. Our 
optimal choice is around 1000 MeV, which is somewhat above the multi-pion production threshold. 
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At this energy, the cross section for quasi-elastic process [v^n — > [i~p or v^p — > fi + n) is 1.5 — 3 
times as large as for the single-pion production cross section. Other channels such as multi-pion 
productions are negligible. 

We note, however, that optimal choice of the neutrino energy may be changed in actual experi- 
ments, depending upon chosen technologies of accelerators as well as detectors. 

The concept of CP-even beam is both detector and neutrino energy dependent, since the and 
cross sections are different for different targets, and the line strengths are different at different 

energies. One may define the CP-evenness 77, using a beam flux T[y^) from ions weighted by the 

cross section o{y^) at definite neutrino energies; 

^; 7) = ^>W-^>fi) (4) 

In order to determine 77 for actual beams, we need to know two quantities: one is the relative in- 
tensity of u e and v e flux, and the other is the relative and cross section. The former quantity, 
J r (i , e )/J-'(u e ), may be estimated from the measured branching ratio of EC and /3-decay. As to the 
latter quantity, cr(z/ /J )/cr(z/ M ), we must include a detector response function (efficiency etc.) in re- 
ality. In this case, the ratio would be best determined in a separate experiment with an identical 
detector component. As an example, let us take 108 Ag 46+ as a beam, and an iron detector. The 
uncertainty of T [y ^ j T [p is ~ 7% at present [2] and CP-evenness r\ « 0.09 selecting the major 
four lines, neglecting energy dependence of cross sections, and including the effect of Fermi motion. 
If the contribution of c(3 integrated over the energy interval ~ 8keV below the threshold is added, 
the CP-evenness 77 « 0. We have so far discussed the idea of a single-ion beam. From the point of 
larger rates, a multi-ion beam is equally interesting. The multi-ion beam consists of a simultaneous 
or time sharing circulation of two different ions of nearly equal Z/A's, each suitable for the bound 
beta and EC, individually having higher rates. For this beam it is important to identify and use 
ions of largest rates for b(3 and EC, respectively. Both candidates of larger bf3 or EC rates than the 
single-ion isotopes are numerous. We only mention two of them; 4g 2 Cd 48+ for b/3, and yjpYb for EC. 
The event rate of neutrino detection depends on the isotope factor of Y/Q with Y the partial b(3 or 
EC decay rate. The isotope 4 22 Cd 48+ , for instance, has the Is b(3 event rate larger by a factor « 12 
than the b/3 rate of 47 8 Ag 46+ , and 7Q 2 Yb EC rate is larger by « 600 than the corresponding EC rate 
of 4? 8 Ag 46+ . The mixed beam made of 4g 2 Cd 48+ (for u e ) and ^ 2 Yb 60+ (for u e ) has similar Z/A, and it 
may be possible to accelerate both of them simultaneously. Their monochromatic neutrino energies 
are 3.0MeV (Is bf3) for 4g 2 Cd 48+ and 5.0MeV (Is EC) for 7o 2 Yb, hence the time sharing circulation 
might be more appropriate. This option of multi-ion beam opens more varieties toward the CP-even 
beam and should be kept in mind for further study. 

After submitting this paper, we have noticed that the measured EC decay rate of hydrogen-like 
i40p r 58+ j on j g about 50% l ar g er than that of helium-like 140 Pr 57 + ion [TT], due to angular momen- 
tum restriction of both leptons and nuclei. [12] [13] A similar mechanism may occur in our system 
i08 Ag 46+ ; iio Ag 46+ ; ii4 Iri 48+ triereDy enhancing our EC rates of Tables 1 and 2. 

In summary, we proposed a new concept to implement the experimental method optimal for de- 
termination of the CP parameter of neutrino mixing by using monochromatic neutrino beams both 
from electron capture and bound beta decay. Feasibility of such beams hinges crucially on production 
and storage of high intensity unstable nuclear beams, and innovative works towards this direction 
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are needed. 
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